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SIGNIFICANCE OF PHOTOGRAPHIC METEOR DATA IN THE DESIGN OF 
METEOROID PROTECTION FOR LARGE SPACE VEHICLES 
by Nestor Clough and Seymour Lieblein 
Lewis Research Center 

SUMMARY 

Described herein are the method of obtaining photographic meteor data in the mass 
-2 0 

range of 10 to 10 gram and the treatment of these data to provide useful engineering 
relations for protecting large space vehicle components against meteoroid impact. The 
assumptions needed to convert these data to meteoroid masses and densities are dis- 
cussed, and the effects of these assumptions on design relations are shown. The inferred 
nature and values of meteoroid density, mass, and velocity are described, and an inter- 
nally consistent meteoroid mass flux model is introduced for the photographic meteoroid 
range. Earth shielding factors as a function of satellite altitude are presented for use 
with the proposed model. Values of average meteoroid density and average meteoroid 
velocity relative to the Earth at the Earth's atmosphere lower than values previously 
used are indicated by an analysis of the available photographic meteor data. Photo- 
graphic meteor astronomy is shown to be a significant source of data in describing the 
meteoroid hazard in space. 


INTRODUCTION 

The space vehicles and power generation systems of the future can be expected to 
have stringent requirements of long duration, large surface areas, and high probability 
of success placed upon them. Components of these flight systems, especially those con- 
taining fluid circuits and life support areas, will have to be protected against critical 
damage due to meteoroid impact. Components with especially large exposed surface 
areas are the propellant tanks of large chemical and nuclear rockets and the waste-heat 
radiators of power generation and coolant systems. 

The critical meteoroid mass, or the mass of the largest meteoroid that the vehicle 
component must be protected against, increases with the vulnerable area, the exposure 



time, and the probability of survival without critical damage (e. g. , ref. 1). This rela- 
tion indicates that large space system components must be protected against meteoroid 
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particles with masses of the order of 10 to 10 gram. Of particular interest to the 
designer of such systems are meteoroid influx rate, density, structure, impact velocity, 
and direction of flight in this particular mass range. Since the design engineer cannot 
predict the exact conditions of these variables, he must rely on the probability of meeting 
specified conditions. 

Numerous satellites and rockets have been launched in an attempt to count directly 

the particle flux in the near-Earth environment (e. g. , ref. 2). The recent experiment, 

Explorer XVI, was the first satellite to measure directly perforation rates in pressurized 

cells. The energy- sensitive detectors of Explorer XVI recorded particle flux rates two 

to three orders of magnitude lower than the rates obtained by previously flown momentum- 
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sensitive experiments for the 10 to 10 gram range (ref. 3). Preliminary reduction of 

these data in terms of a geocentric and a heliocentric distribution of near-Earth particles 

indicates a possible reduction in the apparent discrepancy. 
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Unfortunately, satellite measurements in the 10 to 10 gram range are from four 
to seven orders of magnitude from the range of interest for large space vehicle compo- 
nents, and the data cannot justifiably be extrapolated into this region. An experiment of 
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the Explorer XVI type would be highly desirable in the 10 to 10 gram range. Accord- 
ing to the best currently available estimates of the flux, however, a satellite with approx- 
imately 10^ square feet of vulnerable area and an exposure time of approximately 1 year 

_p 

would yield only one impact of a 10 gram particle. Directly applicable satellite experi- 
ments thus appear to be highly improbable at the moment. 

The design engineer must therefore rely on indirect meteoroid data obtained by as- 
tronomical means (photographic, visual, and radar), which measure particles with 
masses 10” ^ gram and larger. These meteor data fall well within the range of interest 
for the large long- duration space vehicles, and no unjustifiable extrapolation is neces- 
sary. Meteor astronomy is the general term defining the study of meteors. Specifically, 
a meteor is the streak of light caused by a meteoroid traveling at high velocity through 
the Earth's atmosphere. The study of meteor astronomy is therefore specifically limited 
(by definition) to these large particles that leave visible evidence of their passage. Gen- 
eral usage of the term, however, has led to the incorporation of all particles (micro- 
meteoroids to bolides) under this general heading; however, the discussion herein will be 
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limited to photographic -size meteoroids (10 g and larger). 

Considerable meteor data are available in the literature, and various estimates of 
influx rates, densities, and velocities based on these data have been proposed (e. g. , 
refs. 4 to 6). Radical differences in influx rates and densities of meteoroids from pro- 
posal to proposal lead one to be concerned with the significance and accuracy of astro- 
nomically derived predictions. Order of magnitude differences appear in these various 
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proposals, and the design engineer is at a loss to understand the reason for such differ- 
ences and the shortcoming of proposing average density, velocity, and influx rate as uni- 
versally applicable. Use of these various proposals for space vehicle design criteria re- 
sults in large differences in protection requirements and weight among design studies. It 
is important, then, to understand the capabilities and limitations of astronomical reduc- 
tion procedures and how to interpret the data presented in terms of the particle charac- 
teristics and flux. 

One of the objectives of this report is to explain how photographic meteor data have 
been transformed into meteoroid densities and influx rates through the laws of meteor 
physics. The assumptions necessary in the reduction of the data are discussed, and the 
confidence that can be placed on calculated protection requirements is shown. In addi- 
tion, application of the reduced data is discussed and a consistent meteoroid model in the 
photographic region is presented. 


METEOR PHYSICS 
Description 

Meteoritic particles are collectively referred to as meteoroids. Present knowledge 
of meteoroids suggests that they be classified into two major groups, with their origins 
designated as either cometary or asteroidal. The majority of meteoroids (in excess of 
90 percent) are considered to be of cometary origin. These particles represent the de- 
bris left in space after cometary passage. The mechanism of injection and the type of 
particle injected are discussed in a following section. At this point, however, it is well 
to define terms associated with cometary meteoroids. A meteoroid shower is a periodic 
increase of meteoroid flux and represents a group of particles traveling in a definite, 
predictable orbit around the Sun. The fact that many major showers have been corre- 
lated with known comets confirms their cometary origin. Sporadic meteors are not in 
predictable orbits and represent cometary meteors that have had time to depart appre- 
ciably from the initial cometary orbits because of perturbing forces. 

Asteroidal meteoroids are considered to be of different origin from the cometary 
meteoroids. Their origin is believed to be from the asteroids, a number of large 
masses (up to a 800-km diam) in orbit between Mars and Jupiter. As shown in a follow- 
ing section, the asteroidal meteoroids are of a much different structure from the typical 
cometary meteoroids. Their density is about an order of magnitude greater than that of 
the cometary meteoroids, and sometimes they are found to be 90 percent iron. In gen- 
eral, they are a compact stony type of particle, and frequently they survive the passage 
through the Earth's atmosphere and are found on the Earth’s surface. Data are obtained 
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on the larger cometary (and very infrequently on the asteroidal) meteors by photographing 
the meteor trail. 

Ideally, photographic meteor data supply four types of information: particle velocity, 
deceleration, trajectory, and photographic magnitude. These data are obtained by photo- 
graphs of the meteor taken simultaneously with two Baker Super- Schmidt cameras for 
which the base line distance between them is accurately known (approx. 40 km) . In each 
camera is a rotating shutter with a known period, and the area of the sky photographed by 
each camera is known (typically, 5980 sq km at a 90-km altitude) . Unfortunately, not all 
the preceding information can be obtained for all meteors photographed with this equip- 
ment. The fainter meteors and those that do not have a complete trail on at least one of 
the plates cause complications in the determination of the particle characteristics. The 
fainter meteors have very short visible trails, which makes particle deceleration mea- 
surements virtually impossible. 

From the photographic plates of a meteor that is bright enough and the trail of which 
is long enough, the velocity can be calculated by measuring the distance between breaks 
in the meteor trail caused by the rotating shutter. The deceleration is found by measur- 
ing the change in these distances between breaks, and the photographic magnitude is 
found by comparing the meteor brightness with the star background. Finally, each cam- 
era and the line of the meteor determine a plane in space containing the meteor path. 

The intersection of the two planes represents the meteor path in space. Extension of the 
planes allows the determination of the meteor trajectory (i. e. , height and radiant) 

(ref. 7). 


Basic Relations 

An understanding of the physical equations governing the theoretical meteoric pro- 
cess is desirable in order to understand the procedure and assumptions made in calcu- 
lating meteor densities and influx rates from the basic meteor data. Basically, these 
equations include conservation of momentum, conservation of energy, luminosity rela- 
tion, and various combinations of these equations (ref. 8). In order to calculate meteor- 
oid densities and influx rates, however, use of only the conservation of momentum and 
the luminosity relation are required. The equations presented are valid only for solid 
compact bodies traveling through the atmosphere. Fragmentation, or the breakup of the 
meteoroid as it passes through the atmosphere, complicates the use of these equations 
as shown in the following section. 

Conservation of momentum. - Suppose that a compact meteoroid of mass m and 
density p m enters the atmosphere with velocity v. (All symbols are defined in the ap- 
pendix. ) Air particles are trapped on or near the surface of the meteoroid, impart to it 
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the energy of collision, and decrease its forward momentum. The air particles are 
given a forward velocity while the meteoroid suffers a deceleration dv/dt; thus, by the 
principle of conservation of momentum, 

m — = - rs^v 2 (1) 

dt 


where p is the density of the atmosphere, r is the particle drag coefficient, and S is 
the effective cross-sectional area of the meteoroid. One of the variables in equation (1) 
can be eliminated by expressing the effective cross-sectional area S as a function of the 
mass m. An assumption, however, must be made as to the shape of the meteoroid. If, 
for example, the meteoroid is a sphere of density p m , 


S = 





( 2 ) 


1/3 

In general, the meteoroid is an irregularly shaped object, and the term (97t/ 16) ' in 
equation (2) is replaced by a dimensionless shape factor A, which differs from body to 
body. With these transformations, equation (1) can be written as the first fundamental 
equation in meteor theory, called the drag equation: 


dv 

dt 


= -rAp^^m'^pv 2 


(3) 


Lu minosity relation . - The blackbody radiation from the surface of the meteoroid 
represents only a small fraction of the luminosity for the relatively small masses of the 
photographic or visual meteors under consideration. The principal source of the visual 
or photographic radiation occurs in a coma of vaporized meteoric material around and 
behind the meteoroid itself. Practically all the original energy is expended as the cloud 
of meteoric atoms collides with molecules of the atmosphere after the meteoric material 
has left the main body. The luminosity I is therefore assumed to be equal to the energy 
of the mass lost per second multiplied by a luminous efficiency factor r. The luminosity 
is then expressed 


I 


1 dm 2 

_ y 7* 

2 dt 


(4) 


Thus, there are two basic equations, (3) and (4), with mass, velocity, and time as three 
variables; there are four unknown quantities: r, p m , r, and A. Two of the variables, 
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v and t, can be determined from photographic data, and the mass m is still to be found. 
The integration of equation (4) yields 



The terminal mass of the meteor can be shown to be nearly zero, so that 


m = 



( 6 ) 


As shown in a following section, the luminous efficiency r is related to a constant t , 
called the luminosity coefficient, by the relation t = t q v. Then, 



( 7 ) 


The mass referred to in equation (7) is the instantaneous mass at time t (the lower limit 
of the integral) . 

Substituting the mass from equation (7) into equation (3) and rearranging terms 
yields 



Since the mass in equation (3) is an instantaneous value, this equation is valid at a point 
corresponding to the lower limit of the integral. The values of v, dv/dt, and p must 
be the values at this point corresponding to time t. 

With these limitations, the right side of equation (8) can be determined from photo- 
graphic data. The velocity, light intensity, deceleration, duration of the meteor, and 
atmospheric density (determined from the known particle height) are all measured quan- 
tities. In order to determine the density of the meteoroid, t and TA must thus be de- 
termined (or, conversely, if p m and TA are known, the luminous efficiency for a par- 
ticular velocity can be determined). The preceding discussion indicates that meteoroid 
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densities can be calculated readily from equation (8) from photographic data along with 
estimated values of the luminous efficiency and the drag-area coefficient, as long as the 
meteoroids remain compact single bodies. 

Fragmentation . - The typical meteoroid does not meet the conditions of compact body 
flight through the atmosphere. The development of a "correction factor" to be applied to 
the equations of compact meteor flight was therefore necessary. This correction factor, 
called the fragmentation index, was developed by Jacchia (ref. 9) as an answer to the 
"faint meteor anomaly. " This anomaly was discovered when data from the more recent 
Baker Super-Schmidt cameras were compared with the predictions of the "ballistic" 
theory just presented. (The Baker Super- Schmidt cameras are capable of photographing 
much smaller, i. e. , fainter, meteors than cameras previously used.) The durations of 
these meteors were shorter, and their decelerations increased more rapidly than pre- 
dicted by theory. 

These and other observed anomalies can be explained by the process of progressive 
fragmentation proposed by Jacchia. It is theorized that the meteor fragments progres- 
sively along its trail. Small fragments can be detached from the surface of larger me- 
teor bodies without disrupting their unity, but if fragments of similar size are detached 
from small bodies, complete disruption into a cluster of fragments results. This phe- 
nomenon explains the higher frequency of occurrence of anomalies in the fainter meteor 
trails. 

Jacchia developed a quantity, the fragmentation index, as a convenient measure of 
the progressive departure of the observed meteor deceleration from single-body theory. 
The fragmentation index x is empirically determined and is defined as 


X 


d_ 

ds 



(9) 


where a 0 ^ s is the observed acceleration, a^. is the acceleration computed from the 
drag equation (3), and s is the mass loss parameter, defined as 


s = log 




( 10 ) 


where is the original mass of the meteor before it entered the Earth's atmosphere, 

and m is the instantaneous mass of the meteor. Unfortunately, x is not easy to deter- 
mine, because it involves the second time derivative of the velocity; however, x is ap- 
proximately constant for each meteor during the detectable part of its flight (ref. 10). 

The effect of fragmentation on equation (8) can be readily shown by rewriting equa- 
tion (8) as 
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( 11 ) 



This equation was derived, as explained previously, by substituting a value for m into 

the drag equation. The drag equation gives the theoretical value for dv/dt; equations (8) 

/ 2 

and (11) are thus ideal or theoretical equations. In this idealized case, T Q /P m remains 

a constant along the meteor’s path (ref. 10), and, therefore, as long as consistent values 

/ 2 

of v, p, and dv/dt are used, the calculated values of T 0 /P m should remain the same. 
This implies that the observed value of acceleration a Q ^ s is equal to the theoretical 
value of acceleration a^. (as given by the drag equation), and that the fragmentation index 
X as seen from equation (9) is zero. In reality, equation (11) is 



2 

If the value of TApv from equation (3) is substituted, equation (12) becomes 



If m is substituted in the idealized situation from equation (7), 



Equation (14) further shows that if the observed acceleration equals the theoretical accel- 
eration, (t Q / p? a ) equals , which is a constant. 

' 'obs ' ■'ideal 

Differentiating the logarithms of equation (14) yields 



( 15 ) 
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since ( tq / p ^ q ) is a constant for a given meteor, 

yields ideal 


Substituting equation (9) into (15) 



(16) 


and integrating from s Q to s yields 



Equation (17) shows how two values of log ( T 0 /p^ 1 ) on me t eor trajectory correspond- 
ing to two values of the remaining meteoroid mass differ from one another. Again, in the 
absence of fragmentation (y = 0), the two values are the same. 

The calculated value of log (t q /p^j for any meteor can now be related to some ref- 
erence point corresponding to s Q . The data, to be most significant, should be reduced 
to the value of s Q corresponding to the beginning of fragmentation. Unfortunately, this 
value is not generally known, and the data are corrected to a value of s near the begin- 
ning of the meteor. It is indicated intuitively, that the fragmentation index is also an in- 
dex of credibility; that is, the data from meteors with very little fragmentation are 
"corrected" by the least amount and represent data that deviate the least from theoreti- 
cal considerations. 

An excellent example of the use of photographic meteor data with the fragmentation 
correction is given by Verniani (ref. 10). Also discussed is a system of "weights" by 
which the calculated values of log ( t q/p^) at different values of s along the meteor 
trajectory are deemed credible or not. 

With the introduction of the fragmentation correction, it is now generally possible to 
use equation (11) to determine meteoroid density if the luminous efficiency t and the 
drag-area coefficient TA are known. 

Luminous efficiency . - The numerical value of the luminous efficiency t is perhaps 
the largest single influencing factor in the reduction of meteor data. Estimates of t ap- 
pearing in the literature have differed from one another by as much as two orders of mag- 
nitude. Opik (ref. 11) has offered a theoretical evaluation of t for both solid and dust- 
ball photographic meteors. His analysis is lengthy and involved, and by far the most 
complete work to date. The analysis calculates the radiation energy emitted in the wave- 
length range 4500 to 5700 A by iron atoms in a cloud through which a stream of nitrogen 
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Figure 1. - Variation of theoretical luminous efficiency with velocity according to Opik 
(ref. 11). 


atoms is moving with a series of specific high velocities. For the diluted coma condition 
(dust ball), he assumed that the physical process is better represented by single iron 
atoms moving through a cloud of nitrogen atoms, which he shows produces less total ra- 
diation energy. 

The results of Opik's work are shown in figure 1, in which luminous efficiency is 
plotted as a function of meteor velocity relative to the Earth. Note the radical depen- 
dence of t on velocity at values less than 15 kilometers per second. 6pik's work pre- 
dicts a value of zero for luminous efficiency at velocities less than 5. 2 kilometers per 
second. At high velocities, however, the luminous efficiency is very nearly constant for 
both the solid particle and the dust ball. Opik's theoretical values have been used quite 
extensively by investigators in the field of meteor astronomy (e. g. , refs. 12 and 13). 

An attempt has recently been made to measure the luminous efficiency of an artifi- 
cial iron meteor of known mass and density (ref. 14). In this experiment, a seven- stage 
rocket (Trailblazer I) was employed to accelerate a stainless- steel projectile of approxi- 
mately 2. 2 grams into the atmosphere at 10 kilometers per second. The entry of the 
projectile was recorded photographically, and methods of photographic reduction, well 
documented (e. g. , ref. 7), were employed to obtain the luminous efficiency for iron me- 
teors. This value was further reduced to the efficiency of stony particles by assuming 
that the stony particle was 15 percent iron and that all the visible light given off was due 
to iron (ref. 11). This value is then assumed typical of meteoroids that are, for the lack 
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of better information, assumed stony in composition. 

As another source of data, values of luminous efficiency for three photographed as- 
teroidal meteors were calculated by Cook, Jacchia, and McCrosky (ref. 15). Two of the 
meteors were assumed to be stone. The other was identified to be pure iron from its 
available spectrum. Because of the assumed compact nature of the asteroidal particles 
(as opposed to the "fluffy" cometary particles), a good approximation to their density 
was made and the luminous efficiency was calculated with a fair degree of confidence. 

The results of these calculations for the asteroidal particles are shown in figure 2 along 
with the value obtained from the Trailblazer experiment. Fairly good agreement with 
Opik's theoretical calculation for a solid particle is seen. Unfortunately, however, all 
these experimental points are in the low- velocity range of the curve and do not prove, or 
disprove, the constancy of r at high velocities (>20 km/sec) for this type of particle. 

Most investigators have recognized the dependence of the luminous efficiency of a 
cometary meteor upon its velocity. In general, for analytical purposes, an equation of 
the form 

r = T Q v n (18) 

is used to describe this dependence, where t is the luminous efficiency, t q is a con- 
stant (luminosity coefficient), and v is the meteoroid velocity. The exponent n is be- 
lieved to be somewhat less than one (ref. 16). For simplicity in calculation, however, 
and because of the fact that the experimental data can reasonably be fitted within its 
range of error, most investigators use 1 for the exponent n; that is, 

t = t q v (19) 

Shown also on figure 2 is the analytical curve for r according to equation (19) where 

-19 42 -1-34 

t q = 10 ‘ (units of zero magnitude (visual) g cm sec ), given by Whipple as the 

best estimate for use in his calculations (ref. 17). 

Recently, Verniani (ref. 10) made an extensive evaluation of the value of n and of 
the absolute value of t q for cometary meteors. His analysis used 400 precisely re- 
duced meteors and represents the first attempt to analyze existing data on cometary me- 
teors for this purpose. His analysis takes fragmentation into account and is based on a 
knowledge of average meteoroid density. Also, as indicated in the next section, his 
analysis treated so large a sampling of meteors that the effects of unknown individual val- 
ues of the drag- shape factor were reduced. Furthermore, he indicated that n = 1 is a 
good representation of the luminosity law over the complete range of meteor velocities 
(10 to 72 km/sec). It is indicated from the analysis that a linear relation between lumi- 
nous efficiency and velocity (eq. (19)) may be a better representation of the luminosity 
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Luminous efficiency, 



Velocity, v, km/sec 


Figure 2, - Luminous efficiency as function of velocity. Comparison of theoretical, experimental, and analytical 
variations. 
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law for cometary meteors than that proposed by Opik's dust-ball curve. The analytical 
variation for t obtained from Verniani is also shown in figure 2. 

Drag- shape factor TA . - The particle drag coefficient r has been shown to be ap- 
proximately equal to 1 in free molecular flow and approximately 1/2 in continuum flow 
(ref. 16). Levin (ref. 8) has shown that r should assume values of 0. 5 to 1. 1 for condi- 
tions of meteor flight most likely to be encountered. Conditions for typical photographic 
meteors are best represented by r= 1. 1 (free molecular flow) (ref. 10). For the larger 
asteroidal meteors, the conditions of continuum flow prevail with r= 0. 5. 

The dimensionless shape factor A, as defined in equation (2), is a function of the 
shape of the meteoroid. For a sphere, A = 1. 208, and for a cube traveling with one face 
normal to the direction of motion, A = 1. 00. For an oblate spheroid with the minor axis 
in the direction of flight, 



where b is the ratio of the minor to the major axis. Realistic estimates of A based on 
experimental data are impossible to obtain. Probable limiting values of A are difficult 
to determine since a meteoroid can be conceived to have a shape ranging from a long 
rodlike structure to a platelike structure. 

Experimentally, Cook, Jacchia, and McCrosky (ref. 15) found that if TA = 0. 92, the 
average luminosity coefficient of the three asteroidal meteors would be in agreement with 
the Trailblazer (artificial meteor) experiment. The values of the luminous efficiency for 
the three asteroidal particles plotted in figure 2 were calculated by assuming each was a 
sphere (A = 1. 208). An average value for the luminosity coefficient for the three meteors 
was calculated based on the linear dependence of r on v in equation (19), and this value 
was then compared with the luminosity coefficient derived from the Trailblazer experi- 
ment. By adjusting the drag- shape factor of the average asteroidal meteor to 0. 92, exact 
agreement between the derived luminosity coefficients of Trailblazer and the asteroidal 
meteor average was found, and the average value for A was 1. 84. This value implies 
that the "average" shape of the associated asteroidal meteors treated was an oblate 
spheroid with the major axis about 1. 9 times the minor axis. 

In treating a large sample of cometary meteor data, Verniani (ref. 10) indicated that 
the assumption of the meteoroid as a sphere (A = 1. 208) is a good estimate since the 
shape-factor effect will average out if the meteoroid shapes are assumed to be entirely 
random. * Individual calculations of meteoroid density and luminous efficiency are, 
therefore, only as good as the estimated value of the shape factor. For a large sample 

*This assumption may represent a lower limit for A, since a sphere has the least 
ratio of mass to surface area of any solid object. 
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treated as a whole, however, the effects of values of A differing from that of a sphere 
are assumed to ''average*' out, but for individual cases, gross errors in the calculated 
results are possible. 


METEOROID CHARACTERISTICS 
Density 

Until recently, the average density of meteoroids was assumed to be the same as the 

2 

average density of meteorites. This assumption led the astronomer to believe meteor- 
oids were stony iron or iron in nature with a density in the range of 3. 5 to 8. 0 grams per 
cubic centimeter. With this assumption of meteoroid densities, investigators in the field 
of meteor astronomy were able to calculate atmospheric densities from photographic me- 
teor data. With the onset of rocketry came direct measurement of the atmospheric den- 
sity profile at altitudes previously calculated from meteor data. Comparison of the cal- 
culated and measured values of atmospheric density showed that the calculated values 
were consistently low in the altitude range above 7 5 kilometers. This discrepancy of val- 
ues could be resolved if the density of meteoroids were assumed to be much lower than 
that previously used. Whipple (ref. 18), through his icy -conglomerate comet model, pro- 
posed just such a low- density fragile structure for an average meteoroid of cometary ori- 
gin. Opik was also a proponent of the fragile-structure school of thought. His concept 
differs slightly from that of Whipple in that he proposed that the cometary particles enter- 
ing the Earth’s atmosphere and producing meteors are dust balls composed of a great 
number of individual grains (refs. 19 and 20). In effect, both models imply a low overall 
bulk density for a cometary meteoroid. 

Whipple (ref. 13) and others then reversed the process of calculation by assuming at- 
mospheric densities to be known and meteor densities to be unknown. This procedure, 
with improved estimates of the constants involved in the physical equations of meteor phe- 
nomena and recognition of meteor fragmentation, is the one used today. Reevaluation of 
meteoroid densities based on this reversal in calculation led to the proposed low bulk 
density of a meteoroid. 

Whipple, in his latest work (ref. 17), has done this for a group of meteors previously 
reduced by Jacchia (ref. 12). As explained previously, the right side of equation (8) is a 
function only of the photographic data and atmospheric density. Whipple refers to this 

_ O 

value as K . Thus, from equation (8), 

2 A meteoroid is a particle in space; a meteorite is a meteoroid that has survived 
the passage through the Earth's atmosphere and has impacted the Earth. Meteorites are 
generally believed to be asteroidal in origin. 
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( 21 ) 



From the data of reference 12, Whipple obtained an average value of the-right side of equa- 

o ry 

tion (21) based on 88 meteor cases. The value of K = 10 * (units of zero magnitude 
o c 4 _ i q 40 

(visual) g" cm sec ) together with the values of t q = 10 ' and TA = 0. 92 led 

Whipple to conclude that the typical photographic meteor density is 0. 44 gram per cubic 

centimeter. 

There seems to be an inconsistency in the values and therefore in the proposed den- 
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sity of reference 17. The value of r 0 = 10 * is a logarithmic average of the value 

of t determined from the asteroidal particles in reference 15 and the value of r de- 

o o 

termined from the Trailblazer experiment (ref. 14). These values were determined with 
the assumption of TA = 0. 6 (corresponding to a sphere with r = 1/2). As explained pre- 
viously, if TA = 0. 92 is assumed for the asteroidal particles, the values of r Q for the 
Trailblazer and the asteroidal particles could be made to agree. It therefore seems that 
if the logarithmic average of t q is employed, the value of TA = 0. 6 must necessarily 
be used, or if TA = 0. 92 is to be used, the value of t q for the asteroidal particles and 
the value of t q for the Trailblazer experiment must agree with one another and no aver- 
aging is necessary. Apparently Whipple recognized this inconsistency because, on page 29 
of reference 17 he states that a better value of t q is 10 -1 ® - with rA = 0.92. These 
latter values are consistent with the results of reference 15. Employing these latter val- 
ues, Whipple revised his best estimate of meteoroid density to 0. 25 gram per cubic centi- 
meter. 

Verniani (refs. 10 and 21) recently made an exhaustive study of a total of 700 meteor 
cases in order to define meteoroid densities. In reference 10, equation (8) is applied to 
unpublished data obtained by Jacchia. The 359 cases treated in this reference had mea- 
sured decelerations that made the use of equation (8) possible. The meteors were cor- 
rected for fragmentation, and the atmospheric density for application in the right side of 
equation (8) was determined from the known height of the meteor. Since each meteor was 
assumed to be a sphere, TA = 1. 3, and Verniani's derived value of t q = 10"*®' ^ was 
used. 

The results of Verniani's calculations are shown in table I. These average values 
are believed to be significant because of the application of the spherical shape factor, 
which, as explained previously, is assumed to be a valid assumption when applied to a 
large sample. Table I groups the meteors into subsamples of all sporadic meteors, 
sporadic meteors with the aphelion distance Q less than 7 astronomical units, sporadic 
meteors with the aphelion distance greater than 7 astronomical units, and the various 
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TABLE I. - DENSITY OF METEORS COMPUTED BY USE OF 
DRAG EQUATION (REF. 21) FROM UNPUBLISHED DATA 


Meteor 

Number 

Density, 



g/cu cm 

All sporadic meteors 

247 

0.21 

Sporadic meteors with aphelion distance Q < 7 AU 

155 

.23 

Sporadic meteors with aphelion distance Q > 7 AU 

92 

. 16 

Taurids 

23 

.20 

Geminids 

20 

.92 

Aquarids 

16 

.24 

a Capricornids 

13 

. 12 

Quadrantids 

9 

. 15 

Perseids 

9 

.28 

Orionids 

8 

. 14 

77 Cygnids 

4 

. 14 

Lyrids 

3 

.28 

a Hydrids 

3 

.30 

Virginids 

2 

~. 70 

Draconids 

2 

-.001 

Total 

359 



identifiable showers. The separation of the sporadic meteors into these two groups is an 
attempt to show the differences between long- and short-period meteors by defining arbi- 
trarily an aphelion distance of 7 astronomical units as the dividing line. The results 
show that the density is generally smaller than 0. 44 gram per cubic centimeter, with the 
average of all the sporadic meteors slightly greater than 0. 2 gram per cubic centimeter. 
This value agrees with the results as finally derived in reference 17 of an average den- 
sity of 0. 25 gram per cubic centimeter. 

The short-period meteors have a density about 20 percent greater than the long- 
period meteors. Also, the density of shower meteors with small orbits (Geminids, 
Taurids, Quadrantids, and Aquarids) is, on the average, larger than that of showers 
with elongated orbits (Orionids and Perseids). The significance of this difference be- 
tween long- and short-period meteor densities is uncertain. Some investigators believe 
that the short-period meteors have as their origin the inner core of the comet. The 
short-period comet is, on the average, closer to the Sun than the long-period comet and 
will disintegrate more rapidly. The long-period meteors, however, originate from the 
outer surfaces of the long-period comet and may be less compact. 

A search among the meteors studied in this sample revealed no dependence of the 
density on mass or photographic magnitude (ref. 21). Furthermore, the meteors used 
by Whipple (ref. 17) from which he obtained a density of 0. 25 gram per cubic centimeter 
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TABLE H. - DENSITY OF METEORS COMPUTED BY APPROXIMATE 


GROUP METHOD OF CALCULATION (REF. 21) FROM 
PHOTOGRAPHIC DATA OF REFERENCE 24 


Meteor 

Number 

Density, 



g/cu cm 

All sporadic meteors 

280 

0.22 

Sporadic meteors with aphelion distance Q < 7 AU 

203 

.23 

Sporadic meteors with aphelion distance Q > 7 AU 

77 

. 19 

Taurids 

17 

. 17 

Geminids 

16 

.32 

Aquarids 

10 

.49 

Orionids 

12 

. 18 

Perseids 

9 

. 14 

Quadrant ids 

5 

.27 

Bielids 

2 

-.35 

ot Capricornids 

2 

-. 10 

Monocerotids 

1 

-.20 

Draconids 

1 

-.001 

Total 

355 



are of a much greater average mass than those used in reference 21, which produced a 

density of 0. 2 gram per cubic centimeter. This range of mass in which essentially no 

-2 

dependence of density on mass was found is approximately from 100 to 10 grams, or 
four orders of magnitude. 

As pointed out previously, it is not always possible to determine meteoroid decelera- 
tion from photographic plates. This difficulty often makes using equation (8) impossible 
for determining meteoroid density. In reference 21 an approximate method is developed 
for determining meteoroid densities that does not require meteor deceleration but rather 
meteor height at maximum light intensity. The method relates the meteor's shape, 
mass, and density to the atmospheric pressure at maximum light by making certain as- 
sumptions as to the physical strength of the meteor. It is indicated that the method is 
approximate and valid only when applied to a large group of meteors. The method, as 
derived, gives average density directly and loses meaning when applied to single cases. 

For comparison purposes, table n lists values of average density for a group of 355 
meteors (from ref. 21) calculated by this method. The values calculated by this method 
are in close agreement with the direct calculations shown in table I for the sporadic me- 
teors. As is to be expected, the showers calculated by the two different methods do not 
show good agreement with one another when the sample size is small. These results in- 
dicate an average value of meteoroid density typical for sporadic meteors to be about 
0. 2 gram per cubic centimeter. 
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Structure 


As mentioned in the preceding section, meteors, irrespective of size, were earlier 
thought to be compact stones or iron chunks. With better knowledge of the upper atmo- 
sphere and observations of visual and photographic meteors, it was recognized that all 
meteoroids could not have such a high density. Reinterpretation of the photographic me- 
teor data, as explained in the preceding section, has led to the belief that meteoroids 
may have a bulk density of about 0. 2 gram per cubic centimeter. Furthermore, photo- 
graphic data indicate that most meteors break up at the high altitudes at which they first 
appear. Opik (ref. 11) indicates that the visible meteors typically have a crushing 
strength of about 10^ dynes per square centimeter (0. 145 lb/sq in.), or about the crush- 
ing strength of cigar ash. The structure and composition of meteoroids are therefore 
also important considerations. 

The spectra of many meteors have been obtained (ref. 4), but no specific percentage 
compositions are reported. The occurrence of the lines of specific elements in the me- 
teor spectra analyzed, however, is given. The presence of the following elements has 
been detected: sodium, potassium, silicon, aluminum, magnesium, calcium, manga- 
nese, chromium, and iron (ref. 22). Calcium and iron have been recognized as the two 
most prevalent elements in meteor spectra. These data, however, do not allow any con- 
clusions regarding the typical cometary composition because of their qualitative nature. 
In order to obtain quantitative meteor spectra, the luminous efficiency of each element 
in the spectrum must be known. 

McCrosky and Soberman (ref. 14) have indicated another means for estimating the 
composition of typical cometary meteoroids. The relative ratio of heavy elements to 
oxygen in stony meteorites is approximately the same as the relative ratio of heavy ele- 
ments to oxygen in the Sun. Since the meteor spectra indicate the presence of these 
heavier elements, the assumption is made that the same relative ratio of heavy elements 
to oxygen exists in the typical cometary meteoroid as in the Sun and in stony meteorites 
(ref. 14). For this reason it has been concluded that the typical meteor is ’’stony’' in 
composition. 

The observed data indicate then that typically a meteor appears at high altitudes 
(from which low density and low crushing strength is inferred), and the assumption is 
made that its composition is typically stony. The meteoroid structure proposed must 
necessarily produce phenomena in agreement with the observed data and assumptions. 

Two possible explanations exist for the observed phenomena: the first is drag dis- 
tortion and the second is the "dust-ball" concept of a meteoroid. Under the drag- 
distortion interpretation, a compact stony meteoroid entering the atmosphere melts be- 
cause of friction into a spherical drop. Atmospheric drag causes the meteor to flatten 
out into an oblate spheroid with its minor axis in the direction of flight. The increased 
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ablation cross section (frontal area) would then cause the meteor to appear at extremely 
low atmospheric densities. Opik (ref. 11), however, shows that for compact meteors 
under the influence of drag distortion, the length of the meteor trail should decrease with 
increasing brightness. No such effect has been observed statistically, and he therefore 
concludes that the explanation of the increased area of ablation as an effect of deforma- 
tion of a drop is not acceptable. If, then, the compact meteor model deforming under drag 
distortion is discarded, the acceptance of low bulk density, low crushing strength, and 
stony composition as typical of photographic meteors is necessary. It is easy to visual- 
ize this type of particle as being constructed like a cigar ash with layers of soft, loosely 
bound material, or as a stony sponge. 

These cigar ash or stony spongelike particles could have the comets as their origin 
as proposed in the icy- conglomerate comet model of reference 18. In essence, the pro- 
posal is that the comet can be thought of as a large frozen-gas snowball interlaced with 
either the stony sponge or cigar ash material presented previously. Upon repetitive pas- 
sage in proximity to the Sun, the gases boil off and leave the porous, fragile, low bulk 
density material on the outer surface of the comet. Solar radiation, internal pressures, 
and thermal shocks cause the outer surface of the comet to break up into the meteoroids 
observed as meteors. This particular model of meteoroid formation seems to propose 
a particle that satisfies the observed phenomena. 

As mentioned previously, Opik also proposed a meteoroid model to meet the re- 
quirements of low density and low crushing strength (refs. 19 and 20). He proposed a 
model in which the meteoroid is composed of a great number of individual compact grains 
loosely held together; however, he offers no detailed explanation of the injection mecha- 
nisms of such a particle. The differences between Opik's dust-ball meteoroids and 
Whipple’s cometary particles are insignificant, and one could justify the same type of 
origin for the dust ball. In effect, the three meteoroid structures previously discussed 
could be grouped as one: cometary dust balls. 

An interesting proposition has been offered by Donn (ref. 23) for the origin of mete- 
oroids that would give them a somewhat different type of construction. In that reference, 
it is surmised that original meteoroid grains were condensed from the primordial gas of 
the solar nebula when the planets were formed. The grains then agglomerated with a vol- 
atile binding material to form comets. The method of injection of these cometary grains 
into meteoritic particles then is the same as that proposed by the icy-conglomerate 
comet model. 

The significant contributions of the theory of reference 23 is in the proposed struc- 
ture of the cometary meteoroid. It is shown that the particles formed by primordial gas 
condensation would grow as whiskers, filaments, platelets, or some combination result- 
ing in a complex, irregular structure. The original particle can then be described as a 
steel-wool type of structure as one possibility or a many-pointed starlike structure as 
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another. The platelet structure would again propose the ashlike conglomerate. Very 
little change in the original particle structure would result from cometary agglomeration, 
and thus the cometary particles of today would have very nearly preserved the primordial 
condensate structure and composition they originally possessed (ref. 23). This theory of 
meteoroid structure would again fulfill the observed low crushing strength and low den- 
sity typical of photographic meteoroid particles. 

All the previously discussed possible structures for typical cometary meteoroids 
have two characteristics in common: distributed mass with low bulk density and low 
crushing strength. The question then arises as to the impact damage capabilities of this 
type of particle. Existing experimental hypervelocity impact data indicate a very weak 
relation between penetration and particle density (p P^/ (e. g. , ref. 1). These exist- 
ing data, however, were obtained with solid homogeneous projectiles and the density- 
damage relation for heterogeneous, low bulk density type particles is unknown. Further- 
more, the effect of the low crushing strength of the meteoroid on penetration damage is 
also unknown. Again, all experimental and theoretical penetration data to date have been 
obtained with projectiles, the crushing strength of which is orders of magnitude higher 
than the speculated crushing strength of the typical cometary meteoroid. In general, low 
bulk density and low crushing strength in a heterogeneous particle are codependent, and 
thus the individual effects are not likely to be separated. 


Velocity 

A knowledge of meteoroid velocities, besides being required for meteor data reduc- 
tion purposes, is necessary as an input for design protection calculations (ref. 1). 
Ideally, a spatial velocity distribution may exist for meteoroids. To date, a spatial ve- 
locity distribution curve has not been developed. For this reason, ’'average" meteoroid 
velocities have been proposed. It must be remembered, however, that this average ve- 
locity proposal is the velocity relative to the Earth and not necessarily the "impact" ve- 
locity, or velocity relative to the space vehicle in question. For these reasons a discus- 
sion of meteoroid average velocities is felt necessary. 

Figure 3(a) presents a velocity histogram of 285 sporadic meteor cases as reported 
by Hawkins and Southworth (ref. 24). The velocities are relative to the Earth at the 
upper limit of the atmosphere. These data are significant to study since, as discussed 
in a following section, they are the data used by Whipple in proposing his latest influx 
rates (ref. 17) . The distribution of these sporadic data is bimodal and shows a peak at a 
velocity of about 20 kilometers per second and another at about 67 kilometers per second. 
The reduced values between 45 and 55 kilometers per second represent meteors in rela- 
tively unlikely orbits as explained in reference 24. 
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A velocity histogram derived from the data of McCrosky and Posen (ref. 25) is shown 
in figure 3(b). The data consist of 2048 meteor cases that were determined to be sporad- 
ic. The virtual agreement of the two histograms thus verifies the randomness of the data. 
As indicated in reference 25, the method of reduction used for these latter data is not so 
accurate as the method used in reference 24. The randomness of the errors, however, 
validates the treatment of the sample as a group. 

An arithmetic average of the two histograms yields approximately 30 kilometers per 
second for both the data of reference 24 and those of reference 25. This value of average 
meteoroid velocity (relative to the Earth) has frequently been used for the impact velocity 
in meteoroid protection calculations. 

The average velocity obtained from figure 3 represents an average velocity of the ob- 
served sample. The observed velocity distribution, however, is not considered to be 
representative of the true velocity distribution in space because of a velocity bias of ob- 
servation. Such a bias exists because, in general, it is easier to ’’see" the high-velocity 
meteors because they produce brighter trails; that is, relatively more of the existing me- 
teors are observed when their velocities are high. 

As indicated in reference 25, the velocity bias of observation on a photographic plate 
is the combination of two effects. The first effect is seen from equation (4), with 
r = t 0 v. This relation indicates that the light intensity per unit time from a meteor of a 
particular mass is a function of the cube of its velocity. Comparison of light energies 
emitted by individual meteors proportional to the ratio of the cube of their velocities is 
exact, however, only when the masses and duration times are equal (i. e. , when dm/dt 
is constant). 

The second factor contributing to the velocity bias of observation arises from the ef- 
fective exposure time on the photographic plate. The plate exposure time is an inverse 
function of the translational velocity of the point source; therefore, the intensity of light 
recorded on the plate is also inversely proportional to the source velocity. This 
exposure-time effect will influence the observed velocity histograms in a manner oppo- 
site to that of the previous velocity- energy effect, since the faster the meteor travels, 
the less it is exposed on the photographic plate. The velocity histograms are thus 
weighted in the direction of the low-velocity meteors. In reality, however, the situation 
is not quite so simple, since the meteor does not emit light energy as a traveling point 
source. In theory, the visible light energy is emitted by excited meteoric mass parti- 
cles caused by high-energy collisions. These excited particles trail the meteoroid by 
some finite distance and emit as a large number of point sources. The exact dependence 
of exposure time on meteor velocity is further complicated by the existence of a direc- 
tional variation of the velocity vector of the meteor with respect to the photographic plate 
and by the variation in meteor range. (High-velocity meteors have, statistically, a 
greater range. ) 
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The exact observational bias function can be calculated only for individual cases and 
is quite involved. An approximate velocity bias for large samples, however, may be 

O 

chosen as proportional to 1/v by combining the preceding two effects as recommended 
in reference 25. The observed count is then adjusted for this effect. In accepting this re- 
lation for the removal of the velocity bias of the histograms of figure 3, it must be as- 
sumed that, in the photographic mass range of interest (10 to 10 g), there is no pro- 
nounced dependence of meteoroid velocity on mass. The further assumption is made that 
these histograms are representative of any small mass range within the limits of 10 to 
10® gram. Then, for this "constant-mass" histogram there is a higher probability of 
observing the high-velocity particles than the low-velocity particles. 

In establishing a more nearly correct count, it is assumed that all existing particles 
are observed for velocities between 70 and 75 kilometers per second in figure 3, and a 
correction for the velocity bias of the photographic observations proportional to 1/v is 
made for the other velocity ranges. On this basis, the "true" distribution of particle 
velocities for the samples considered assume the shapes shown in figure 4. As expected, 
these histograms indicate many "unobserved" particles of low velocity, and significantly 
this indicates that there is more flux present than observed. The average velocity of the 
true distribution of the data of reference 24 (fig. 4(a)) is 20. 2 kilometers per second and 
of the data of reference 25 (fig. 4(b)) is 20. 3 kilometers per second, values that are sig- 
nificantly less than the nominal 30 kilometers per second of the "observed" distributions. 


MASS INFLUX RATE 

The design engineer must know how many meteoroids of a particular mass will be 
encountered by the vehicle during its trajectory. The determination of meteoroid influx 
rates from photographic meteor data has been thoroughly treated by Hawkins and Upton 
(ref. 26). The classical approach employed in reference 26 is generally accepted and 
will be briefly explained herein. The photographic data obtained, as explained previ- 
ously, yield the light intensity and photographic magnitude of the meteor directly. An 
examination of these data shows an increase in the number of the faint meteors as com- 
pared with the bright meteors, which indicates an increase in the number of meteors 
with decreasing mass. 

Hawkins and Upton (ref. 26), from a random sample of 285 sporadic meteor trails, 
derived an equation for the cumulative influx rate N, expressed as a function of the 
photographic magnitude M over the range of photographic magnitudes from -2. 0 to 

r 

+4. 0, in the form 


where the constant a is the slope of the curve of log N against and N Q is the num- 
ber of meteors with magnitudes equal to or less than zero. Here, the photographic mag- 
nitude Mp is an average value over the meteor trail. By simply counting the number of 
meteors with magnitudes less than or equal to zero, log Nq was determined: 

log N q = -4. 34 (23) 

-2 - 1 

(where Nq is in km hr ). This number was based on an average rate of observation 
by the viewing cameras of 0. 278 meteor per hour over a collection area of 5980 square 
kilometers of sky. The number of meteors increased by a factor of 3. 4 per magnitude 
and thus produced a value of a of 0. 537. Equation (22) then becomes 

log N = 0. 537 M - 4. 34 (24) 

Ir 

The representation of meteor flux rate in terms of photographic magnitude given by 
equation (24) is perhaps the most accurate form of this type of information. It is uncon- 
taminated by further calculations and is strictly a derived representation of raw data. 

The accuracy of the equation is a function only of the accuracy of measurement of the 
photographic magnitude and the area of the sky photographed. This representation of me- 
teoroid influx rate, however, is not very usable for space-vehicle-protection require- 
ments, since a knowledge of meteoroid mass-influx rates is necessary. 


Hawkins and Upton Rate 

Hawkins and Upton (ref. 26) indicate that they also calculated individual masses for 
each sporadic meteor in their sample of 285 sporadic meteors. This calculation can be 
made, as explained previously, by the use of equation (7) with an assumed value of the 
luminosity coefficient. The choice of a particular value for the luminosity coefficient 
then implies that a consistent value of an average meteoroid density is determined. This 
value is approximately 3 grams per cubic centimeter for the influx rate of Hawkins and 
Upton. They found that the cumulative influx rate N could be expressed as a function of 
the meteoroid mass m in the form 

log N = aj log m - log Nj (25) 

where the constant a^ is the slope of the curve of log N against log m and Nj is the 
number of meteors with masses equal to or greater than 1 gram. Determination of the 
constants a j and log N^ yielded 
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log N = -1. 34 log m - 2.64 


(26) 


p 1 

(where N is in km hr ). Equation (26) is a more useful form of the meteoroid influx 
representation for a space-vehicle designer. It is not so accurate a representation of the 
data as equation (24), since it does have implicit in it an assumed luminosity coefficient. 


Whipple 1963A Rate 


Equations (25) and (26) are the basic equations used by Whipple (ref. 17) in determin- 
ing his proposed mass- influx rate, referred to as the ’’Whipple 1963A flux rate. " The 
basic difference between this proposed influx and equation (26) as used by Hawkins and 
Upton is the luminosity coefficient used in calculating meteoroid masses. Since this co- 
efficient is the basic difference, Whipple adopted equation (26) and corrected it to a new 
luminous efficiency in the following manner: Equation (7), with r = t q v, shows the mass 
m as 



(7) 


Since the masses of Hawkins and Upton (ref. 26) differ from the mass of Whipple (ref. 17) 
only by a different luminosity coefficient, the ratio between the two mass values is 


m W _ T o, H 
m H T o, W 


or 


m w = 



(27) 


where the subscripts W and H 


refer to Whipple and Hawkins, respectively. 


Then 


T 0 TJ 

log m w = log — ^ — + log m fl 
T o,W 

Substituting log m H from equation (28) into equation (26) yields 


( 28 ) 
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log N = 1. 34 log 


(29) 


'Y 

- 1. 34 log m - 2. 64 

T o,W 

which can be used to correct the Hawkins and Upton equation to any value of the luminos- 
ity coefficient. 

Hawkins and Upton used r Q = lO -20, 91 . Whipple (ref. 17) proposed r Q = 10“ 19 ' 42 ; 
equation (29) therefore becomes, for the Whipple values, 

log N = -1. 34 log m - 4. 63 (30a) 


(where N 


-2 - 1 

is in km hr ) or 


log N = -1. 34 log m - 14. 18 (30b) 

(where N is in m sec ). Equations (30a) and (30b) (differing only by their units), give 
the proposed influx rate of Whipple (ref. 17), the Whipple 1963A flux rate. The relations 
of equations (30) do not contain the Earth shielding factor of 1/2 that Whipple included in 
his proposed influx rate. The reduction in flux by a factor of 1/2 as proposed by Whipple 
(ref. 17) is not realistic for the general case, since obviously the effect of Earth shield- 
ing is a function of the altitude of the vehicle, and can be 1/2 only when the vehicle is at 
the effective meteoroid boundary of the Earth T s atmosphere. 

Included in the proposed influx rate of equations (30) is an average meteoroid den- 
sity consistent with the assumed luminosity coefficient. This value is 0. 44 gram per 
cubic centimeter. The influx rate of equations (30), with an average velocity of 30 kilo- 
meters per second and an average density of 0. 44 gram per cubic centimeter, has be- 
come the most widely used meteor flux model to date. 


Revised Rate 

As discussed in three previous sections, better values of meteoroid mean density 
(based on a new luminosity coefficient) and meteoroid average velocity appear to be avail- 
able. It seems reasonable, therefore, to propose a meteoroid flux model based on these 
values. In general, equation (26) can be corrected to a new value of the luminosity coef- 
ficient t q by using 

in -20. 91 

log N = -1.34 log m - 2.64 + 1.34 log — (31) 
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Verniani (ref. 10) arrived at t q = 10 and P m = 0. 2 gram per cubic centi- 

meter as representative of cometary meteors. If this value of r Q is accepted as the 
best estimate to date, equation (31) is corrected to yield 

log N = -1. 34 log m - 4. 23 (32) 

(where N is in km hr ). As shown previously, the average velocity of the observed 
sample (fig. 3, p. 21) is about 30 kilometers per second. The use of equation (32) with 
an average velocity of 30 kilometers per second and an average density of 0. 2 gram per 
cubic centimeter would yield a consistent flux model for the observed sample. 

As discussed previously, the velocity bias of observation indicates that many parti- 
cles are unobserved. Since the design engineer is more interested in the actual flux dis- 
tribution than in the observed distribution, equation (32) should be further corrected to 
account for these inferred "unobserved" particles. The mass influx rate derived by 
Hawkins and Upton (ref. 26) has in it a correction of observation based on the assumption 
that all the large mass particles are observed. The following observational bias removal 
assumes that all the high-velocity meteors were observed, and that regardless of mass 
the low-velocity meteors went largely unobserved. 

The velocity- corrected distribution yielded v c = 20 kilometers per second compared 

with the observed distribution value of v ^ of 30 kilometers per second. The decrease 

in average velocity was attributed to the inferred increase in the number of low-velocity 

particles. This increase in particle number was assumed previously to be proportional 
o 

to 1/v for constant particle mass. The increase in sample number, for essentially 
constant mass, would thus be proportional to the square of the ratio of the sample aver- 
age velocities; that is, 


or 



N obs \ v c / 


A log N = 0. 352 


(33) 


Equation (32) is then modified to yield 

log N= -1.34 log m - 3.88 (34) 

- 2 - 1 

(where N is in km hr ). Equation (34) with p m = 0. 2 gram per cubic centimeter and 
v = 20 kilometers per second is then a consistent meteoroid flux model for a near- Earth 
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Figure 5. - Deduced influx rates based on data of Hawkins and Upton (ref. 26). 


environment representing an estimate of the "actual" rather than the observed distribu- 

-2 0 

tion. The mass range over which this model can be assumed to apply is 10 to 10 gram, 
which is the range of the observed data. 

The various flux models discussed herein along with their associated values of den- 
sity and velocity are shown in figure 5. The flux models of Hawkins and Upton, Whipple 
1963A, and the solid curve of equation (32) are all representative of the observed refer- 
ence data, and the differences are wholly attributed to the assumed luminosity coefficient. 
The dashed curve is the inferred "actual" distribution as discussed previously in equa- 
tion (34). 

The revised actual distribution (dashed curve) indicates an increase in the meteor 
hazard as compared with the Whipple 1963A curve. The proposed associated lower den- 
sity and velocity, however, counteract much of this increase, as shown in a later section. 


EARTH SHIELDING FACTOR 

When a flux model is used for meteoroid protection calculations for an Earth satel- 
lite, the influence of the Earth as a shield should be considered. Whipple recognized 
this influence but applied a single factor of 1/2. In reality, the Earth shielding factor E 
is a variable and can be so computed as a function of the altitude of the satellite above 
the Earth T s surface that, at any altitude, 

N = EN (35) 

where N is obtained from equation (31). 

In calculating E, it will be assumed that the meteoroid population is isotropic. Fig- 
ure 6 represents a unit hemisphere of space of unit radius, with its center at the satel- 
lite; r is the true radius of the Earth, h is the satellite altitude above the surface of the 
Earth, and R is the "effective" radius of the Earth. Using the effective radius of the 
Earth recognizes that most meteors burn out in the upper atmosphere at altitudes of 70 
to 100 kilometers. The atmosphere of the Earth therefore provides an effective shield 
to meteoroids and thereby serves to increase the radius of the occluding body. For 
practical purposes, the atmosphere at 90 kilometers, or approximately 50 nautical miles, 
should provide an effective shield to meteoroid passage; that is, Ar = 50 miles. With the 
use of this assumption, the effective Earth shielding as a function of altitude can be de- 
rived from reference to the geometry of figure 6. 

The shielding factor of a satellite at altitude h above the Earth's surface is deter- 
mined to be 


30 






( 36 ) 


E _ 1 , V(r + h) 2 - (r + Ar) 2 

2 2(r + h) 

Figure 7 shows the Earth shielding factor as a function of altitude based on the as- 
sumption that the atmosphere is an effective shield up to 50 nautical miles. The Earth’s 
occlusion is of importance to satellites in orbit at less than about 2500 nautical miles. 


DESIGN IMPLICATIONS 


From the preceding discussion it is obvious that the reduction of photographic meteor 
data into meaningful information useful to the space-vehicle designer is not a simple and 
precise task. Furthermore, there are many limitations and assumptions involved in the 
reduction methods that can have notable effects and implications on the design protection 
requirements. These limitations and assumptions can best be appreciated by examining 
a typical relation used for predicting required armor thickness to prevent critical mete- 
oroid damage. 

Reference 1 presents an equation for required single-sheet armor thickness as a 
function of mission parameters and meteoroid hazard as 
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In terms of only meteoroid properties, 


fcr 


(38) 


where p m 
factors a 


and v are the meteor bulk density and impact velocity, respectively, 
and -]3 are related to the meteoroid flux according to the equation 
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The exponent -/3 is the slope of the general curve of log N against log m (eq. (25)) and 
a is the same as N^. 


32 



Meteoroid Density 


In a previous section the calculation of the density for a particular meteoroid was 
shown to be dependent on several factors other than the data obtained from the photo- 
graphic plate. In particular, the luminosity coefficient r Q and the drag-area coefficient 
TA, as well as the correct application of the fragmentation index, were shown to be highly 
influencing factors. The uncertainty introduced in the value of the meteoroid density by 
t he uncert ainty in these factors is easily seen from equation (11) as proportional to 

Vt 0 CTA) 3 . It is impossible, however, to place uncertainty limits on either t q or TA, 
since at best they are educated estimates. The treatment of large groups of data, as ex- 
plained previously, is assumed to nullify the uncertainty in the drag-area coefficient and 
thus to reduce the uncertainty in the "average" meteoroid density primarily to a function 
of t q (if fragmentation is assumed to be correctly applied). Fortunately, equation (38) 
shows that the required armor thickness for space-vehicle protection bears a very weak 
relation to the meteoroid density. This weak relation, which was derived for solid homo- 
geneous particles, may be a conservative relation for the low-bulk-density, low-crushing- 
strength meteoroid particle. 


Meteoroid Velocity 

The treatment of estimated impact velocity warrants further study. The attempt 
herein to remove the observational bias of the data has resulted in a reduction of the 
average velocity of the meteor sample relative to the Earth from a nominal 30 to 20 kilo- 
meters per second. The correction presented is not precise but rather is an attempt to 
show an effect that can markedly influence the total meteoroid hazard. If an average me- 
teoroid velocity relative to the Earth is to be better defined, a more precise treatment is 
needed in which the observational bias is removed from individual meteor cases by con- 
sidering all the factors involved and the fact that all meteors are not of the same mass. 
The treatment presented assumed that the velocity histograms were independent of mass. 
It is probable that large particles will be observed regardless of their velocity; this pos- 
sibility makes the correction presented invalid for this region of mass. 

The value of 20 kilometers per second deduced herein as a typical or average mete- 
oroid velocity is relative to the Earth at the limit of its upper atmosphere and is not nec- 
essarily the value that should be used in equation (38). Ideally, the value in equation (38) 
should be the effective meteoroid impact velocity relative to the space vehicle for its 
mission. Uncertainties in relative impact velocity affect the required single- sheet armor 
thickness drastically (proportional to v 3 ^ 3 ), and a significant effort is therefore indi- 
cated for the determination of accurate values. Meteoroid relative impact velocity is a 
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function of the vehicle trajectory in space and the spatial meteoroid velocity distribution. 
The second of these two items is at present undefined. Much work remains in the defini- 
tion of this spatial velocity distribution. The heliocentric as well as the geocentric rela- 
tions involved are necessary, and a method of incorporation of the results into protection 
equations such as equation (38) must be devised for a precise analysis of the hazard. To 
date, meteor data are lacking for the definition of this distribution. For an Earth- orbiting 
vehicle at high altitude (long periods), a 20-kilometer-per- second relative impact velocity 
for particles in the photographic mass range may be an adequate estimate. More advanced 
missions, however, such as those that leave an Earth orbit or that involve orbits with 
short periods or high eccentricity, will require a more sophisticated analysis of the esti- 
mation of meteoroid impact velocities. 


Mass-Influx Rate 

The influence of the mass-influx rate on the required armor thickness in equation (38) 
is indicated by the values of a and /3. 

The value of 0, the slope of the curve of log N against log m, equal to 1. 34 estab- 
lished by Hawkins and Upton in 1958 (ref. 26) has remained unchanged. Although a more 
precise value of j3 can probably be determined from an increased sample mass range 
and sample number, j3 = 1. 34 can be assumed to be sufficiently valid for use over the 
photographic mass range. 

Uncertainties in a will reflect uncertainties in 6 approximately proportional to the 
1/4 power (if /3 = 1. 34 is assumed). The value of a is seen from equation (31) to be di- 
rectly related to the uncertainty in the luminosity coefficient t q raised to the 1. 34 power. 
Comparing the influx rates of Whipple (eq. (30)) and Hawkins and Upton (eq. (26)) shows 
a two-order-of-magnitude difference in the value of a (from 2. 64 to 4. 63) and corre- 
spondingly in the value of flux. This variation is wholly attributable to the differences in 
the values used for r . This indicated relation reflects the strong dependence of the re- 
quired armor thickness on the absolute value of r . As explained previously, this abso- 
lute value is difficult to determine; therefore, confidence in t q (and finally a) becomes 
related to the "educated guess" made in determining r . 

In addition to uncertainties in a (related to the accuracy of t q ), the value of the flux 
will depend on the magnitude of the velocity bias of the observation and the effective re- 
moval of this bias from the observed sample. The dashed curve of figure 5 (p. 29) indi- 
cates an "actual” distribution inferred from observational data based on an approximate 
removal of the observational bias, as explained previously. This representation of an 
"actual" rather than the observed influx rate is currently felt to be more meaningful for 
protection calculations. A better curve could be inferred from a larger sample of data 
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by removal of the bias from individual meteor cases. The "probability" of observation, 
therefore, does not have to be made independent of mass, as explained previously. 

Finally, the influx rates discussed herein should specifically be limited in use to a 
near- Earth environment. The effect of gravity focusing by the Earth is implicit in the 
measured flux rate and therefore in the inferred actual flux rate. Since the focusing ef- 
fect tends to concentrate meteoroid particles in the immediate vicinity of the Earth, flux 
rates could be significantly different from the observed values at large distances from the 
Earth. The focusing effect, being velocity sensitive, requires the treatment of individual 
meteor cases for greater precision in the determination of the flux. 


Comparison with Whipple 1963A Hazard 


The meteoroid hazard as proposed herein based on the inferred "actual" flux 
and revised values of density and velocity can be compared with the Whipple 1963A hazard 
on the basis of equation (38). The ratio of the required armor thickness for the revised 
hazard 5^ based on the dashed- curve environment model proposed herein to the armor 
thickness given by the Whipple 1963A hazard model 6^ for the same mission is given by 



(40) 


The following table lists the various input values for each model: 


Model 

Meteoroid 

Average 

Slope of mass 

log a. 


density, 

meteoroid 

against influx 



Pm’ 

velocity, 

rate, 



g/cu cm 

V, 

0 




km/sec 



Whipple 196 3 A 

0.44 

30 

1. 34 

-4.63 

Revised environment 

.2 

20 

1. 34 

-3.88 

(dashed curve, fig. 5) 






Substituting into equation (40) the various values of p , v, /3, and a for each model 
yields 


—5. = l. 03 
5 W 


(41) 
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The revised hazard indicates that about 3 percent more armor is needed than predicted by 
the Whipple 1963A hazard model. This minor increase in required armor thickness is 
due to the counterbalancing effects of an increased indicated flux and a decreased average 
impact velocity. For all practical purposes the two models yield the same armor re- 
quirements. 


Meteoroid Structure and Shape 

Two major factors have arisen in the determination of the average meteoroid density 
that can have a further marked influence on the required protection. These factors re- 
late to the shape and the structure of the meteoroids and the basis of available impact re- 
lations. Current knowledge of penetration and damage effects due to the impact of hyper- 
velocity particles (e. g. , eq. (38)) has been obtained from experimental and theoretical 
studies based on solid homogeneous particles, usually spherical or cylindrical (equal 
length and diameter) in shape, and with relatively large crushing strengths. Current es- 
timates of the meteoroid properties, as described herein, however, indicate that mete- 
oroid particles may be substantially different from the man-made particles in both struc- 
ture and shape. 

First, it is believed that meteoroids may vary considerably in shape and range possi- 
bly from long needlelike rods to flat plates. Experimentally, the shape of a solid impact- 
ing particle and its orientation at impact have been shown to have a major effect on the 
resulting penetration (e. g. , ref. 27). If these same effects persist for meteoroids, sig- 
nificant variations in impact damage may result from individual meteoroid particles of 
the same nominal characteristics, that is, mass, density, and velocity. 

Secondly, and perhaps more important, meteoroid particles, by virtue of their known 
constituents, low bulk density, and possible origins, are believed to possess relatively 
low crushing strength with a structure containing a heterogeneous distribution of mass. 

If this is indeed true, then the resulting impact effects could be different from the cur- 
rently available solid- particle relations, and, hopefully, less severe. Furthermore, if 
such differences exist, it would be important to determine the possible influence on these 
differences of such factors as meteoroid velocity, shape, and mass. Theoretical and ex- 
perimental investigations may shed some light on possible impact differences, while ad- 
ditional meteor astronomy techniques or spacecraft experiments may be useful for fur- 
ther structure and shape determinations. 


CONCLUDING REMARKS 


It is shown herein how photographic meteor data and methods of analysis assembled 
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from various scientific sources are evaluated and interpreted for use in the design and 
development of meteoroid protection for large spacecraft components. It was indicated, 
however, that considerable uncertainties remain in astronomically derived meteoroid 
properties. Since the pertinent meteoroid properties for the larger particles, which con- 
stitute a significant threat to space devices with large vulnerable areas, cannot be obtained 
by practical flight experiments, further effort should be directed toward the reduction of 
the uncertainties in the photographic meteor analysis process. Efforts are needed to in- 
vestigate the impact characteristics of meteoroid-type particles as discussed herein and 
to define more accurately meteoroid structure, shape, and composition. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 7, 1965. 
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APPENDIX - SYMBOLS 


A 

\ 

a 

a l 

a obs 

a t 

b 

C 

D 

D 1 

E 

h 

I 


K 


m p 


m 

m 

m 


e 

oo 


N 

N c 

N 

I 

N„ 


obs 


shape factor 

vehicle vulnerable area 

slope of magnitude influx rate 

slope of mass influx rate 

observed meteor acceleration 

theoretical meteor acceleration 

ratio of minor axis of meteoroid 
particle to major axis 

sonic velocity of impacted target 

constant in penetration equation 

constant 

Earth shielding factor 

altitude of satellite above 
Earth’s surface 

intensity of light given off by 
meteor 



-1/3 

_dv) /J. dt 
L \pv 2 dt / J v ? 

photographic magnitude 
meteoroid mass 

meteoroid mass at end of meteor 

meteoroid mass before it enters 
Earth’s atmosphere 

cumulative influx rate 

corrected cumulative influx rate 

observed cumulative influx rate 

Earth- shielded cumulative in- 
flux rate 


Nq cumulative influx rate of zero- 

magnitude meteors 

Nj cumulative influx rate of 1-gram 

particles (same as a) 

n exponent 

P penetration 

P(0) probability of no critical damage 

Q aphelion distance 

R effective radius of Earth 

r mean radius of Earth 

Ar difference between effective ra- 

dius and mean radius of Earth, 
R - r 

S frontal area of meteor 

s mass loss parameter 

s Q mass loss parameter at particu- 

lar point on meteor trail 

T mission time 

t time 

t time at end of meteor 

e 

v meteoroid velocity 

v average meteoroid velocity 

v c corrected average velocity 

v 0 bg observed average velocity 

a cumulative influx rate of 1-gram 

particles (same as Nj) 

/3 slope of mass agsinst influx rate 

(1.34) 

r particle drag coefficient 
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6 

e 


p 

P m 
Pi 


required armor thickness to pre- t luminous efficiency 

vent critical damage Tq luminosity coefficient 

constant in penetration equation fragmentation index 

( 2 / 3 ) 

Subscripts: 

atmospheric density 

H Hawkins 

meteoroid density 

R revised 

target material density 

W Whipple 
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